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, Abstract 

, We compute the Green's function for the Hodge Laplacian on the symmetric spaces M x E, 

where M is a simply connected n- dimensional Riemannian or Lorentzian manifold of constant 
curvature and E is a simply connected Riemannian surface of constant curvature. Our approach 
, is based on a generalization to the case of differential forms of the method of spherical means and 

on the use of Riesz distributions on manifolds. The radial part of the Green's function is governed 
by a fourth order analogue of the Heun equation. 
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Our purpose in this paper is to compute the Green's function for the Hodge Laplacian in some special 

classes of symmetric spaces of Riemannian or Lorentzian signature, more precisely those spaces which 

are products M x S, where M is a simply connected n-dimensional Riemannian or Lorentzian manifold 

of constant curvature and E is a simply connected Riemannian surface of constant curvature k [16J. 

[ The spaces we consider are thus not of constant curvature and they are not conformally flat, unless 

one of the factors reduces to a point. They contain in particular the Robinson-Bertotti solutions of 

the Einstein-Maxwell equations, as well as the de Sitter and anti-de Sitter spaces, the latter spaces 

' corresponding to cases in which S reduces to a point and M is of Lorentzian signature. 

' 

' The approach we take is based on a natural extension to the case of differential forms of the classical 

method of spherical means [12] , which has been applied with great success to the computation of the 
Green's function for the scalar wave equation in symmetric spaces [1^. The method of spherical 
means has the advantage of leading naturally to the use of Riesz distributions, which are fundamental 
analytical tools in the study of linear wave equations on manifolds. Thus we avoid the need for special 
ansatze on the form of the action of the Hodge Laplacian on the Green's function, as is often the 
case in the Physics literature [2 [11] . In the case of spaces of constant curvature, Riesz distributions 
provide a rigorous way of reducing the computation of the Green's function to the study of an ordinary 
differential equation governing its radial part. We will see that a similar reduction occurs in the case 
of the product spaces considered in this paper. 

The computation of the Green's function for the Hodge Laplacian in de Sitter and anti-de Sitter 
space is a problem of considerable interest in Physics, where differential forms appear as basic variables 
in various supersymmetric field theories, including supergravity. In this context, the Green's function 
corresponds to the propagator of the theory, which is the basic object of interest. To our knowledge. 
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the earliest results on this problem appear in papers by Allen and Jacobson |2j and Folacci [TUj ■ Their 
results have been used in the recent works of a number of physicists, e.g. [51 ll3[ [5]. The case of product 
manifolds seems to have received less attention in the literature and presents challenges of its own, 
as we shall see below. Beyond the interest of this problem in Physics, our motivation has been to 
illustrate the power of the method of spherical means as applied to operators acting on sections of 
vector bundles, and to show that it can be adapted efficiently to the case of product manifolds. 

One of the main features of our result is that in the case in which neither E nor A/ is reduced 
to a point, the expression for the Green's function of the Hodge Laplacian acting on p-forms (with 
l^p^n— 1) has its radial part governed by a fourth order ordinary differential equation with four 
regular singular points, that is a "higher order Heun equation" , which appears to be new. This is in 
contrast with the case in which S reduces to a point, where the ordinary differential equations for 
the radial part of the Green's function (which are coupled by the Levi-Civita connection terms) can 
nevertheless be decoupled by a suitably chosen differential substitution into second-order equations 
of hypergeometric type. In all cases, we obtain a closed form expression for the Green's function as a 
bundle-valued distribution. 

There are some natural extensions that one could pursue for the results that we have obtained 
in our paper. Of immediate interest would be the extension of the method of computation of the 
Green's function for the Hodge Laplacian to warped products of spaces of constant curvature. One 
could also develop the formalism of our paper in the case of the Laplacian on spinors, along the lines 
of [3] either in the plain or warped product cases. A more challenging objective, which would perhaps 
be of less interest in Physics, would be to adapt to the case of differential forms of arbitrary degree the 
analysis carried out by Lu |15| for the Green's function for the Laplacian on 1-forms over symmetric 
spaces which are not maximally symmetric, such as the complex hyperbolic spaces and some classical 
symmetric domains of higher rank. 

Our paper is organized as follows. In Section[21 we briefly recall the definition and key properties of 
the spherical means operators on differential forms defined on Riemannian space forms, as well as their 
relation to the action of the Hodge Laplacian. Section [3] contains a concise account of the method of 
spherical means in the Lorentzian case, and introduces the representation of the fundamental solution 
of the wave equation in Lorentzian space forms in terms of Riesz distributions. In Section 01 we 
consider the case of Riemannian products, where we establish in Theorem 14.51 the expression of the 
Green's function for the Hodge Laplacian in terms of the spherical means operators and the solutions 
of a fourth-order ordinary differential equation with four regular singular points, that is a "fourth- 
order Heun equation" which has not been studied in the literature as far as we know. An important 
element in our proof is a spectral decomposition result (Proposition |421) for the second-order operator 
of singular Sturm-Liouville type corresponding to the radial part of the Hodge operator on surfaces, 
which we prove in Section [71 In Section [51 we establish using Riesz distributions the expression of the 
Green's function for the Hodge Laplacian in the case of Lorentzian products. This is done in Theorem 
15.31 and also involves a fourth-order radial equation of Heun type which arises now from the use of 
Riesz potentials. Finally, in Section [51 we brieffy show for the sake of completeness that in the case in 
which the surface S reduces to a point, that is when the manifold is a simply connected Riemannian 
or Lorentzian space form, the Green's function obtained through our approach is expressible in terms 
of the solutions of a second-order equation of hypergeometric type, as shown in ^ and fTl^. 

2 Spherical means on differential forms defined on Rieman- 
nian space forms 

In this section we shall briefly review some results on spherical means for differential forms defined on 
the simply connected Riemannian space form M of constant sectional curvature k. Proofs and further 
details can be found in Ref. [1^. We denote by BM{x,r) the ball centered at the point x G M oi 
radius r and set S]\i{x,r) — 9i?Af(a;,r). It is well known that 6*^/(2;, r) is diffeomorphic to a sphere 
for any < ?■ < diam(Ai^) and that these spheres foliate M\{x} (minus the antipodal point of M is 
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positively curved). We use the notation dS* for the induced area measure on S'M(a;,r), whose total 
area is ^ 

m{r):=\S{x,r)\ = \S-'\(^^^^ . (1) 

Throughout this paper we shall use the determination of the square root with nonnegative imaginary 
part, which has a branch cut on the positive real axis and is holomorphic in C\[0, +00). We also 
recall that the injectivity radius of M coincides with its diameter, which is +00 if fc ^ and -k j^J —k 
if /c < 0. The space of smooth p- forms (resp. of compact support) on M will be denoted by ^^(M) 
(resp. %{M)). 

Let p e C°(M X M) be the distance function on M. For any degree -p we define double differential 
p- forms Tp,fp by setting 

To(a;, a;') := 1 , t\{x,x') := ddl p[x,x') , (2a) 

Vk 

Tp{x,x') := -Tp-i{x,x') AA' Ti{x,x') for 2 < p ^ n (2b) 

and 

fo(x, x') := , Ti(x, x') := —dp{x, x) A' p{x, x) , (2c) 

Tp{x,x') := Ti{x,x') AA' Tp-i{x,x') for 2<p^n (2d) 

at the points {x, x') E M x M where the distance function is smooth. In the above formulas accented 
and unaccented operators act on x' and x, respectively. Globally, this defines bundle-valued distribu- 
tions Tp,Tp e P'(M X M,AP{T*M) m AP(T*M)), where Ap{T*M) K1 Ap(T*M) is the vector bundle 
over M X M whose fiber at (x, x') is Kp{T*M) ® Kp{T*,M). 

Given a smooth section * of Ap{T*M) M Kp{T*M) and uj e fiP(M), we introduce the notation 

^{x,x') ■uj{x) := *{^{x,x') A*(j{x)) , 
^{x,x') ■' oj{x') := *'{^{x,x') A' *'uj{x')) 

for the left and right inner products r(AP(T*M) K Ap(T*M)) x Op(M) ^ C°°{M) «) nP{M). Left 
and right products can be also defined for ^/ gV(M x M, Ap{T*M) M Ap{T*M)) and w € nP{M). 

An important property of the double differential forms defined above is that if x' does not belong 
to the cut locus of x, then the sum Tp{x, x') + fp(x, x') provides the parallel transport operator for p- 
forms along the unique minimal geodesic connecting x and x' . More precisely, consider luq G Ap(T*M) 
and let 7 be the minimal geodesic with 7(0) = a; and 7(1) =x'. Then [Tp{x,x') + Tp{x,x')]-wo = a)(l), 
where the smooth map ui : [0, 1] Ap{T*M) satisfies the differential equation for parallel transport 
along 7, which with a slight abuse of notation can be written as 

V^(t)'^W = , w(0) = Wo • 
In particular, Tp{x, x) + fp(.T, x) is the identity map in Ap{T*M). 

Definition 2.1. Let < r < diam(M). The (Riemannian) spherical means of a smooth j3-form 
u) e QP{M) on the sphere of radius r are defined as 

Mrujix) := / Tp(x, x') ■' w(x') dS(x') , (3a) 

m'T) Js{x,r) 

MrU}(x) / fp(x, a;') •' lo{x') dS{x') . (3b) 

m{r) Js{x,r) 
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Example 2.2. If a; G n^{M), then A4r(^ ~ and Air^^ coincides with the usual spherical mean on 
constant curvature spaces [T4j , which is customarily written as 



MMx) = J^^_^ ivix + re) dsie) 

in Euclidean space. If a; e Q"{M), Mr^ = and *{M.ri^) = 

We shall hereafter denote by Am dd* + d*d the (positive) Hodge Laplacian in M and consider 
the measure on (0, diam(M)) given by d^M{i') '■— m{r) dr. From the point of view of their applications 
to the Hodge Laplacian, the key property of the spherical means is given by the relation [T^ p. 47] 



A 



diam(Af) 



/(r) Mri^ + f{r) Mri^j d^Af (r) 

diam(A/) 



(^F{f, f){r) MrUJ + M^uj) dy.M{r) , (4) 



where uj S r2^'(M) and /, / are arbitrary smooth functions for which the latter integrals converge. The 
functions f),F{f, /) are given by 



F{fJ){r) LMf{r) + pk {2 cac^ Vkr + n - p - l)/(r) - 2cosVIrcsc^ Vkr f{r) 
F{f, f) [r) := LMf{r) + {n - p)fc [(2 csc^ Vkr +p- l)/(r) - 2 cos Vkr csc^ Vkr f{r 
Lm being the radial Laplacian 



q2 g 
Lm TTT + (n- l)Vk cot Vkr — 



(5a) 
(5b) 

(6) 



This formula remains valid for / S L\^^{{Q^ diam(M)), d/is) if ^ is compactly supported and the above 
derivatives and integrals are understood in the sense of distributions. 



3 The method of spherical means in the Lorentzian case and 
Riesz potentials 

In this section we shall review the extension of the method of spherical means to Lorentzian spaces of 
constant curvature and some useful properties of Riesz potentials. Again we refer to [12] for details. 
We shall denote by M the simply connected Lorentzian manifold of constant sectional curvature k and 
signature (— , +, • • • , +). The Minkowski (fc = 0) and de Sitter spaces (fc > 0) are globally hyperbolic, 
but the anti-de Sitter space (fc < 0) is not (it is, however, strongly causal |16j). For this reason it is 
convenient to consider a domain D d M which will be fixed during the ensuing discussion. When 
k ^ 0, D can be taken to be the whole space M, whereas when fc < we shall assume that D is a 
geodesically normal domain of M [12 , that is, a domain which is a normal neighborhood of each of 
its points. 

We introduce the notation J^{x) for the points in D causally connected with x, i.e., the points 
x' £ D such that there exists a future directed, non-spacelike curve from x to x' . Furthermore, we set 

5+(x,r) :={.T'e J+(x):p(x,x')=r} 

and call dS its induced area measure. Obviously 

JU^)=[J S+{x,r), with 5+(a;,r)n5+(a;,r')=0 if r^^r'. (7) 

We also use the notation p : D x D ^ [0, +oo) for the Lorentzian distance function. 
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Definition 3.1. The Riesz potential at x with parameter a is the distribution ^ e 'D'{D) defined 
for Re a ^5 n by the integral 

Rl ^ a / p(x, xT'-^cfix ) Ax' , (8) 

where (j) is an arbitrary smootli function compactly supported in D, dx' stands for the Lorentzian 
volume element and we have set 

For the sake of completeness, let us recall that the map a i— > Rf) ^ is holomorphic in C provided 
that a 1-^ Rf) ^{(j)) is an entire function for any smooth function with compact support in D. A 
basic result in the theory of Riesz potentials is the following theorem [171 112] ■ 



Theorem 3.2. For any x Cz D, the map a i— > Rf) ^ can be holomorphically extended to the whole 
complex plane. Moreover, Aji/i?^ ^ = ^d^x '^'^'^ x ~ ^x- 

Corollary 3.3. For any 4> G d^{D), u{x) :— Rj^ ^{cf)) solves the hyperbolic equation Amu — (j) in D. 

We define the bundle- valued distributions Tp,fp over D using the Lorentzian distance p as in 
Section [21 If a; e _D and x' £ J^i^) does not lie on the causal cut locus of x, Tp{x,x') + ■fp{x,x') is 
again the parallel transport operator Ap{T*D) Ap{T*,D) along the unique future-directed minimal 
geodesic from x to x' . In particular, Tp(x, x) + Tp{x, x) is the identity map on Ap{T* M). 

Definition 3.4. Let uj £ nf^{D). Its (Lorentzian) spherical means of radius r are defined as 



M^uj{x) -^-^ / r(x, x') ■' w(a;') dS{x') , 



m{r) 



Sl(x.r) 



M^oj{x) := .y;;^ / t{x, x') •' Lj{x') dS{x') . 

S+{x,r) 



m{r) 

where x £ D and m is given by Eq. ([T]). 

We shall henceforth drop the scripts D when there is no risk of confusion. It should be noticed that 
the latter operators do not exactly describe a spherical mean in the strict sense since in the Lorentzian 
case m(r) does not yield the area of S'^i^x, r), which can be infinite. However, the fundamental relation 

Aa/ j (/(r) MrUJ + /(r) MrOj) dfiM{r) = J f)(r) MrUJ + F(f, f){r) Mrto) dpM{r) , 

(10) 

also holds true in the Lorentzian case, with F(/, /) and F{f, /) defined by Eq. (O. 



4 Riemannian products 

In this section we shall solve the Poisson equation 

AAfxsV' = ^, (11) 

oj being a compactly supported differential form in the Riemannian product space M x S. If Af x S is 
compact, we assume that lo is orthogonal to the harmonic forms of the same degree in M x S, which 
is the necessary and sufhcient solvability condition of the latter equation jTS] . 
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It obviously suffices to analyze Eq. pTjl when u £ ilP{M) (g) r2^(E) for some integers p and q. 
If we denote by Qm and f^s tlie Riemannian volume forms in M and S, respectively, the above 
orthogonality conditions amount to imposing 



a;f]MAOs = if (p, g) - (0, 0) , ^ w A = if (p, g) = (0, 2) , (12a) 

w A l^s = if (p, (?) (n, 0) , / w = if (p, g) = (n, 2) (12b) 



when M X E is compact. 

We begin by discussing these exceptional cases. The explicit solution of the scalar equation (p = 
q — 0) for functions w in M x S (and more general symmetric spaces) is well known and can be 
found in [T3]. Concerning the other exceptional cases, let us denote by and the Hodge star 
operator of M and E and observe that *j:uj, or *^{*m^) are C°° scalar functions on M x E, 

respectively, when {p,q) = (0,2), (n, 0) or (n, 2). As *m and commute with the Hodge Laplacian 
of M X E and preserve the orthogonality relations (fT2|) . for (p, q) e {(0, 2), (n, 0), (n, 2)} Eq. (fTTj) turns 
out to be equivalent to the scalar case by Hodge duality. Therefore we can henceforth assume that 
(Pi q) ^ {(0, 0), (0, 2), (n, 0), (n, 2)}. For concreteness we shall also assume that k and k are nonzero; 
the case where one of the curvatures is zero can be treated along the same lines (and is considerably 
easier). 

We find it convenient to introduce the notation Ss and Sg for the spherical means of radius s in 
the Riemannian surface E and call 

, , ^ sinv^s , 
dME(s) ■■=2Tr ^ds 

the radial measure on (0, diam(E)). It should be observed that the decomposition 

n*{M X E) = np{M) (g) r2«(E) 

defines a natural action of the Laplacians and spherical means operators of M and E on fl*{M x E). 

A simple but important observation is the following 
Lemma 4.1. Let (3 e l^^(E). Then 

/'diam(E) /'diam(S) 

As/ w{s){Ss+Ss)l3dfij:{s)^ L,«;(s) (5, + 5,)/3d/is(s) , 



where w € i?iQj.((0, diam(E)), d/is) anc 



92 g 

Lo = £2 -TTT - V«cot V«s J (13a) 
OS'' OS 



LMs) ■■= Low{s) + ^^^""^^^ . (13b) 

1 + cos WKS 



Proof. We saw in Example 12.21 that Sgu; — when q — 0. Hence in this case we easily obtain the 
desired formula from Eq. ([4]) by noticing that Low — F{w,w). When q = 1, this follows directly from 
Eq. ([4]) since Liw = F{w, w) = F{w, w). The case q = 2 is essentially equivalent to g = by Hodge 
duality and uses that L2W — F{w, w). □ 

The formal differential operators Lq given by define self-adjoint operators (which we still 
denote by Lq for simplicity of notation) with domains given by the functions u G H^{{0, diam(E)), d/is) 
such that LqU £ i^((0, diam(E)), d^s), Hm su'(s) = 0, and 

slO 

lim ( diam(E) — s) u'{s) =0 if k > and q = . 

s^diam(S) 
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Next we state a useful proposition on the spectral decomposition of the operators Lq. The proof, which 
consists of an application of the Weyl-Kodaira theorem and some manipulations of hypergeometric 
functions, is postponed until in Section [71 

Proposition 4.2. Let c^^q :— ~jK if k < and c^^.q := 2KSqi if k > 0. Then for each g = 0,1,2 
there exist a Borel measure Pq on [c^^q, oo) and a (/is x pq) -measurable function Wq : (0,diam(I])) x 
[cK,g,oo) —* such that: 

(i) Wq{-,X) is an analytic formal eigenfunction of Lq with eigenvalue X for pq-almost every A £ 

[CK,g,00) 

(a) The map 

/'diam(E) 







U{s) Wq{s, •) dps{s) 



defines a unitary transformation Uq : L^((0, diam(I])), d/is) ^ -^^([cn.g, oo), dpg) with inverse 
given by 

Uq^U:^ / Wq{;X)u{X)dpq{X). 

(Hi) Lf g{Lq) is any bounded function of Lq, then 

g[Lq)u= / g{X)Wq{-,X)UqU{X)Apq{X) (14) 

in the sense of norm convergence. 

Remark 4.3. Explicit formulas for all the functions involved are given in Lemmas 17.11 and 17.21 By 
definition, pq is supported on the spectrum of the self-adjoint operator Lq. 

We shall construct a solution to Poisson's equation pT|) as 

= j Wq{s,X){Ss+Ss){f{r,X)Mr + f{r,X)Mr)oo<lp^{s)ApM{r)dpq{X) , (15) 

where the integral ranges over [ck,^, oo) x (0, diam(M)) x (0, diam(I])). By construction, LqWq{s, X) — 
Xu}q{s,X) for pg-almost every A, so from Lemma 14.11 and the fact that Aa/xsV' = ^m'>P + '^sV' it 
follows that 

Aa/xsV-^ / Wq{S + S)[{F{f,f) + Xf)M + {F{f,f) + Xf)M\iudpj:dpMdpq, (16) 



with F and F given by ([5]). 

Let us consider the system of coupled ODEs 

F{fj) + Xf = 0, ]?(/,/) + A/ = 0, (17) 

depending on the nonnegative parameter A. Let us introduce the variable z :— sin^ which takes 
values in (0, 1) when r e (0, diam(Af)) and M is compact and in (— cx3,0) when M is noncompact. 
With an abuse of notation we shall temporarily write f{z) or f{z) for the expression of the functions 
/ or / in this variable, that is, for /(-^ arcsiny^). This allows to write Eqs. ^Tl} as 



z{l-z)r{z)+'-^{l-2z)f'{z) 
z(l-z)/"(z) + ^(l-2z)/'(z) 



^+(n-p)(^^+p-l 



2z[l — z) 

(18b) 
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with (. := X/k. We can solve (|18ap for / and substitute in the second equation to arrive at a hnear 
fourth order differential equation for /, namely 



(19) 



with the rational functions Qj being given by 
n + 4 n + 4 4 



Q2{z) 
Qi{z) 



; - 1 z 1 - 2z ' 

16(n^ + 2{p + 3)n - 2/ + + 3))z'* + 32(2/ - ~ 2{p + 3)n 
2{£ + 3))z^ + 8(3^2 + 5(p + 3)n + 3n - Sp^ + £ + 4(£ + 3) + 6)2^ + 8(?j2 - 
ip + 3)n-3n-e- 6)z + 71^ + 612 + 8 j (l6z^ - 48z^ + h2z'^ - 2^z^ + 4^^ 

2pn^ + - 2p^n + 4pn + 2ln + 2n - Ap^ + 4£ + 4( - 2prr' - 



2pn - 2ln - 2n + 2p^ + 2p-2()z + 4(2pn^ + n' 
2p + 21) z"^ I (^Az^ - lOz'' + 8z^ - 2z2) , 

2np''' + n^p^ - 2ip^ + p^ - 2np + 2nlp - 2p + f 
2n£p -f)z + 4(/ - 2np^ + n^p 



2p^n 



2p^n + 2pn + 2£n + 2n - 2p^ 



2V 



21- 
2 



4{-p- 
2£j3^ -p^ + 2nlp 



2np^ 



2„2 



n p 



4z^ - I2z^ 



13z*-6z-' + z' 



The fourth order equation pop plays a crucial role in the computation of the Green's function of 
the Hodge Laplacian. A simple computation shows that it has four regular singular points at 0, i, 1 
and 00, so it can be understood as a fourth order analogue of the Heun equation. For our purpose it 
is important to make the following observation. 

Proposition 4.4. For all £ ^ there exists a unique real solution $+(z,£) of Eq. (fT9)) with the 
asymptotic behavior at 



<^+{z,£) 



4z^ 



k{n-2)\S'' 



— +0{z''-^) ifn^S and $+(z,^) 



47r 



+ 0(1) if n^2, (20) 



and the fastest possible decay as z tends to ~oo. If £ > or p ^ 0,n and £ ^ 0, there is also a unique 
real solution $_(z,£) of the latter equation such that z"5"~^$„ (z, £) is analytic in the closed interval 
[0, 1] and has the asymptotic behavior (|20p at 0. 

Proof. Eq. (|19p is a fourth order Fuchsian equation with poles at 0, ^, 1 and 00. Suppose that n ^ 3. 
Then characteristic exponents at and 1 are —^,l — j,0 and 1, so that must be given by 



$+(z,£) = 



fc(n-2)|S"-i| ^^-^^"^ + cMz) + c,Mz) 



where (j)^ stands for the local solution of p9)) which is analytic on (— oo,0) and asymptotic to z'^ at 
0. The characteristic exponents of at 00 are ^[n + 1 ± {{n + 1 — 2p)^ — M)^/^] and i[n — 1 ± ((n — 1 — 
2p)^ — 4£)^/^], so that each local solution ipi, has well-defined asymptotics at 00 and the real constants 
Ci and C2 can be chosen so as to obtain the fastest possible decay. When n = 2, the local solutions at 
zero behave as z~^, log |z|, 1 and z and the same reasoning applies mutatis mutandis. 

The analysis on the interval (0, 1) is similar. The characteristic exponents at | are 0, 1, 3 and 4, so 
all the solutions of the equation are analytic at this point. The function <i>_ must therefore be given 

by 

4 



<i>_(z,£) 



fc(n-2)|S" 



— (j)l^ll{z) + Ci(t)o{z) + C2(/)l(z) , 
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where the local solutions (j)^, are now analytic in (0, 1) and asymptotic to at 0. The constants ci 
and C2 should now be chosen so as to ensure that <&_ does not have any terms asymptotic to or 
z^~T . With some effort it can be shown that this is always possible when p 0,n or £ ^ 0, which is 
not surprising since these conditions are directly related to the solvability of the associated equation 
{k£ + Am)^' = a, with M compact, a € QP{M) and M^Q. □ 

Now we have all the ingredients to prove the main result of this section. 

Theorem 4.5. Let us set f{r,X) :— Wq(0, A) <i>signfe(sin^ A/fc) and define a function f{r,\) by 
means of Eq. (jl8ap . i.e., as 

f, LMf{r,\)+pk{2csc^Vkr + n-p-l)f{r,X) 

f(r, A) :— ^ j= -= . 

2pk cos V kr csc^ v kr 

Then the function ip given by solves Eq. ^ for uj G QP{M) ® r2«(S). 

Proof. We start by noticing that, by Proposition l4.4l and the restriction on the possible values of (p, q), 
the function (f>sign fc (sin^ A/fc) is well defined for all values of r and A in the integration range. It 
should be noticed that the definition of /(•, A) ensures that it has the same asymptotic behavior at 
that /(•, A), namely ^ 



(n-2)|S"-i|r"-2 

An important observation is that the distribution defined by A), /(•, A)) d/is (which equals 

F{f{-, A), /(•, A)) d^s) is in fact Wq{0, A) times the Dirac delta supported at 0, since 

^diam(M) ^ 

f{r,X)LM'p{r) +pk(^{2csc^Vkr + n-p-l)f{r, A) 

-2cosVfcrcsc2 VA?r/(r,A)V(r)jd/is(r) = Wq{0,X)ip{0) (21) 



for all (p e C^([0, diam(Af)]). This immediately stems from the fact that F{f{-, A), /(•, A))(r) is zero 
for all r by construction (cf. Proposition 14. 4p and the asymptotic behavior of / and / at zero. In 
particular, as A4o + A4o is the identity map it easily follows that 

f.diam(M) 

[{F{f, f) + Xf)Mr + (F(/, /) + Xf)Mr]^d^lM{r) = . 



As a consequence of this, Eq. reduces to 



M = J wiO,X)w{s,X){Ss+Ss)iudfiT.{s)dpg{X). (22) 

When the integral in A ranges over [c^,^, oo), it immediately follows from Proposition l4.2l and Lemma [7.3l 
on the pointwise convergence of the integral (fT4|) that 



AV' = (5o +So)uj = u, 

and the claim in the statement follows. □ 

Remark 4.6. Being constructed using only the distance function and the metric on M and E, the 
symmetric Green's function that we have constructed is equivariant under the isometrics of M x S. 
Remark 4.7. Some comments on the uniqueness of the solution to Eq. pT|) are in order. When M x S is 
compact, it follows from standard Hodge theory that for the above values of {p, q) we have constructed 
the only solution to the equation. The case when M or E is noncompact can be analyzed using 
that ker(AMxs) = ker(AM) ® ker(Ax;) and that the kernel of the Hodge Laplacian of the hyperbolic 
m-space acting on s-forms is {0} if n 7^ 2s and infinite dimensional otherwise [7l[8]. Hence for the 
above values of {p,q) Theorem 14.51 also yields the only solution to (fTTjl when fc > 0, k < and 
(Pil) ^ {(0, 1), (fi, 1)}, when fc < 0, k > and {p,q) ^ {(f , 0), (|, 2)}, and when fc < 0, k < and 

(p,g)y^(t,i)- 
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5 Lorentzian products 



In this section we shall solve the equation 



%jj = UJ 



(23) 



for a compactly supported form lo e QP{D) (g) Vt'^{Yj) by constructing an advanced Green's operator 
for Amxs- We assume that M and E respectively have Lorentzian and Riemannian signature. As 
in Section [31 if M is not globally hyperbolic [k < 0) we restrict ourselves to a geodesically normal 
domain Z? x E. If M is globally hyperbolic {k ^ 0) we simply set D := M . The symbols Mr and Mr 
will stand for the Lorentzian spherical means in D, but other than that we will use the same notation 
as in the previous section. We recall that for any globally hyperbolic manifold Z) x E there exists a 
unique advanced Green's operator ^^{D x E) — > f7*(_D x E) for the Hodge Laplacian j^. 

As discussed in the previous section, when (p, q) e {(0, 0), (0, 2), (n, 0), [n, 2)} Eq. ((23)) is equivalent 
by Hodge duality to the scalar wave equation, whose Green's function is well known [T3]. Hence we 
shall assume that (p, q) does not take any of the above values and that both k and k are nonzero. 

The results in the previous section and the definition of the Riesz potentials strongly suggest that 
we analyze the system of ordinary differential equations 



A), /„(•, A)) (r) + A /„(r, A) = C„_2 r"-"-^ 
^'(/a(-, A), /„(., A)) (r) + A /„(r, A) = C0.-2 r""""' 



(24a) 
(24b) 



where A is positive, Ca is given by Q and we assume for the moment that Re a > n + 2. We 
find it convenient to write this equation in the variable z := sin^ writing fa{z) or fa{z) for the 
expression of the functions fa{r,X) or /Q(r, A) in terms of z with some abuse of notation. Defining 
the fimction ha{z) : 



siny^)" " ^, Eqs. (pi)) now read 



z(l-z)/;'(z) + -(l-2z)/;(z) 



z(l-z)/;'(z) + -(l-2z)/;(z) 



2z(l - z] 



+ n — p — 1 



faiz) 



p{l - 2z) ^ 



£+{n~p) 



2z(l - z) 



p-1 



2z{l - z)' 
Uz) 



(25a) 



h.{.)- ^''-'?}'-''h .{z), (25b) 



where £ A/fc. We can combine Eqs. 



2z(l-z) 

5]) to obtain a single fourth order equation for f^, namely 



dz* 



3=0 



(26) 



where the rational functions Qj were defined in Section 31 

Ha{z) := qo{z)ha{z) + qi{z)h'^{z) + q2{z)h'^{z) 

and 

(z - l)z (4nz2 + 8z2 - 4nz - 8z + n + 4) 



92(2) 
90(2) 



2(z-l)^z^ 
P(2z-1) ' 



p(2z- 1)2 



2z {£{z- 1)(1 - 2z)2 + p(-zp + p + n(z - 1) + z)(l - 2zf - 2z + 2) 



p{2z-lf 

Using the results from the previous section it is not difficult to prove the following 
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Proposition 5.1. Under the above hypotheses, Eq. (|26p has a unique solution which is 

continuous in [0, 1] if k < 0, and a unique solution ^a.~{z,£) which is continuous in (— oo,0] and has 
the fastest possible decay at infinity ifk<0. These solutions are real for real a. 



Proof. Let us assume that n ^ 3 and fc < 0. We omit the discussion of the case n ~ 2, which goes 
along the same hnes by replacing z^"^ by log \z\ in the discussion. We saw in Proposition 23] that the 
homogeneous equation (|19p has four regular singularities at 0, ^, 1 and oo. As Ha is continuous on 
(— oo, 0] for a > n + 2, the method of the variation of constants yields a particular solution $o of pH)) 
which is analytic in (— oo,0) with possible singularities at of order and z^~t . The function in 
the statement of the theorem is thus given by 

- + Ci(f>^^{z) + C2(l)l^^{z) + C^Mz) + CiMz) (27) 

in the notation of Proposition 14.41 The constants ci and C2 are chosen so that 

limz5"^$a = 0, 

zTO 

while C3 and C4 are chosen so as to obtain the fastest possible decay at infinity. 

Suppose now that fc > 0. Let us observe that for a > n + 2 the function Ha is continuous in 
[0, i) U (i, 1], whereas it diverges as (z — i)"^ at ^. However, as the characteristic exponents of the 
homogeneous equation at ^ are (0, 1, 3, 4), it is standard that the method of the variation of constants 
yields a particular solution $0 of Eq. (|26p which is continuous (actually, analytic) at ^ . The desired 
solution is obtained, in the notation of Proposition 14.41 as 

'^a,+ iz,i) = ^oiz) + Ci0_.j(z) +€201-^(2) + C30o(-z) + C40l(2;) . 

Here the constants Cj are chosen so that 

limz5"^$Q.+(z,£) = lim(l - z)* "^$a.+(z, £) = 0, 

zjO zj.0 

i.e., so as to remove the singularities of order ^ and ^ — 1 at and at 1. From the proof of Proposi- 
tion it stems that this determines the constants cj. □ 



The AP(r^M)-valued distribution defined for each x e D as 
^",d(0) ■=Ca /r"-"(X,+A?,)0(x)dAiAf(r) 



Ca pix, x' {t{x, x' ) +t{x,x')) ■' (t)ix')dx' , ^el^g(M), 



is easily seen to be holomorphic for Re a > n, with fixed x and A. By Theorem 13.21 the function 
a '—^ TZ" u admits an entire extension, and as t(x,x) + t{x,x') = idAp(T*A/) it follows from this 
theorem that 

<d(</')=0(^)- (28) 
In what follows we shall set fa{r,\) — <&a,signfc(sin^ -^^,A/fc), with $q,± as in Proposition 15.11 



and define fa so that Eq. ([24]) holds true, that is, as 

iAf/a(r, A) + pk{2 csc2 ^/kr + n-p- l)/„(r. A) - C„_2r° 



/a(r,A) 



2pk cos ^/kr csc^ \/kr 
It is natural to define another vector-valued distribution by 



K:d{4>) ■■= / (/a(r, A) Mr + fcAr, A)M,)^(a;) d^^A^r) , G n^oiM) 
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for Re a > n — 2. By Eqs. PH)) . and (|28p it follows that the two distributions that we have just 
defined are related by 

In fact, from the above equation and Proposition (STT] it is standard [IZlll] that the function a ^ 'R-x'd 
can be holomorphically extended to the whole complex plane. Thus we are led to a useful generalization 
of Corollarv l3.3l that the summarize in the following 

Proposition 5.2. The distribution 7?."'^ is a holomorphic function o/ a G C for fixed x ^ D, 

A G spec(Lq). Moreover, for any cj) G ^Iq{D) the differential form ^{x) :— T^"'jj{<P) solves the equation 
{Am + X)^ = (j) m D. 

The action of 7?."'^ naturally defines a bundle-valued distribution on M x S, which we do not 
distinguish notationally, which acts on the first factor of any oj e flP{D) (g) r2'?(S) to yield an element 
of AP(T*M) (g) r2«(S). Namely, if uj{x, y) = uji{x) ® ^2(2/) this action reads as 

'K%{^){y) ■= Kloi^i) ® ^2(2/) , 

with (x, y) E D X I]. This result allows us to express the solution of Eq. as follows. 
Theorem 5.3. The differential form 



ip{x,y) := J Wg{0, X) Wq{s, X)n"ljj{uj){y) dfiM{r) dnj:{s) dpg{X) , 

where the integral ranges over [cq^^, 00) x (0, diam(S)) x (0, diam(A/)), solves Eq. (I23p for any compactly 
supported u G nP{D) (g) 



Proof. If immediately follows from the definition of Wq and Proposition 15.21 that 



AMxY.ip{x,y) = I Wq{0,X)wq{s,X) {X + AM)'R-"'D{Lo){y) d^M{r) dfi^{s) dpq{X) 



Wq{0, A) Wq{s, A) uj{x, y) d^s(s) dpq{X) . 

As in the proof of Theorem 14.51 now Proposition 14.21 and Lemma 17.31 show that 

AAfxsV'la;, y) = {Sq + Sq)uj{x, y) = uj{x, y) , 
completing the proof of the statement. □ 

6 Constant curvature spaces 



We saw in Eqs. p7|) and ((24)) that the radial behavior of the Green's operator in the product manifolds 
M X S is controlled by a fourth order analogue of the Hcun equation. On the other hand, it is well 
known [21 [11] that the Hodge Green's function on M (as well as the scalar Green's function of any 
simply connected rank 1 symmetric space |14j ) is controlled by simple hypergeometric equations: this 
merely reflects that the geometry of M x S is considerably more involved than that of its individual 
factors. For the sake of completeness, in this short section we shall discuss the simplifications that 
give rise to hypergeometric functions when one considers the equation on Af, that is, when the factor 
E collapses to a point. 
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We shall therefore analyze the following system of coupled ODEs: 



z(l-z)/"(z) + -(l-2z)/'(z) 



z(l-z)/"(z) + -(l-2z)/'(z)- 



1 



2z(l - z) 



+ n — p — 1 



1 



2z(l-z) 



p(l - 2z) 

p-l) />) = 

(n - 2z) 



h{z) 



2z{l - z) 



(29a) 



/(z) . (29b) 



When h{z) — this equation controls the radial part of the Green's operator for the Poisson equation 
on M X S studied in Section |4l whereas when h is the function ha defined in Section [5] we obtain the 
equation for the radial part of the retarded Green's function. It is natural to look for functions a, d, b 
and b (possibly depending on p and £) such that the functions 

g{z) := f'iz) + pa{z)f[z) - pb{z)f{z) , g[z) f'{z) + (n - p)d{z)f{z) - {n - p)b{z)f{z) (30) 

satisfy a first-order system of ODEs (equivalent to (HH)) of the form 

z(l - z)g'[z) + {n- p)A{z)g{z) + pB{z)g{z) = h{z) , (31a) 

z(l - z)g{z)+pA{z)g{z) + {n - p)B{z)g{z) = h{z) . (31b) 

When this reduction can be performed, it is possible to express the functions / and / in terms of 
the solutions to the decoupled second-order equations for /, / and g, g that one obtains from ([50]) 
and (PT|) . As before, by Hodge duality and the fact that the scalar case is well known allows us to 
assume that p =/= 0,n. 

After some manipulations one finds that these equations amount to imposing the following condi- 
tions: 

{z ~l)zb{z) - B{z) = 0, 
{z ~l)zb{z) - B{z) ==0, 
n(i - z) + {n - p){z - l)za{z) - pA{z) = , 
- z) + p{z - l)za{z) + {p - n)A{z) = , 



1 



1 



1 



2^^__^+-)+(p- n)A{z)b{z) + {p- n)B{z)a{z) + {z - l)z b'{z) = , 

+ -] - pA(z)b(z) ~ pa(z)B(z) + (z- l)zb'(z) = , 
2 \z — 1 z J 

1 I 

— — + 1 - n + p H h (p - n)a{z)A{z) + {p - n)B{z)b{z) + (z - l)z a'{z) = , 

2(z — l)z p 



1 



l-p 



- pa{z)A{z) - pb{z)]3{z) + (z - l)za'(z) = . 



2z(z — 1) " 71 — p 

Thus we have eight equations with eight unknowns, but solving the above system is in general a 
formidable task. 

The crucial observation is that this system of ODE simplifies considerably when ^ = 0. In fact, in 
this case we can set a = d and b = b, which immediately leads to the solution 

1 — 2z - 1 

a(z) = a(z) = — ■ , biz) = biz) = . 

^ ' ^ ' 2z(l - z) ' ^ ' ^ ' 2z(l - z) 

Thus one derives that when £ = the equations satisfied by the new functions g and g are 

z(l - z)g\z) + {n-p){^- z)g{z) + ^pg{z) = h{z) , (32a) 

z(l - z)g'{z)+p{l - z)g{z) + \{n - p)g{z) = h{z) , (32b) 
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and by isolating g in (|32ap one readily finds that 



z(l - z)g"{z) + (I + 1) (1 - 2z)g'{z) - {n - p){p + l)g{z) = H{z) , 

with H{z) := h'{z) + Shi^. This is a hypergeometric equation, which can be solved in terms of 
associated Legendre functions by the variation of constants method. 

Now it suffices to notice that the equations controlling the radial behavior of the Green's operators 
on M are obtained from those of M x S by collapsing E to a point, which amounts to setting £ = 0. 
Thus from this discussion and the previous sections we recover the (non-rigorous) result of AUen- 
Jacobson and Folacci [TT] that the Green's function for the Hodge Laplacian in a simply connected 
space of constant curvature can be expressed in terms of hypergeometric functions, as in the scalar 
case. Detailed albeit somewhat formal discussions can be found in the aforementioned references; 
details are omitted. It should be stressed that for arbitrary values of £ the radial equation (|29aP do 
not seem to admit an analogous reduction. 

7 The radial equation on surfaces 

In this section we shall prove Proposition 14.21 and derive explicit formulas for the elements appearing 
in the statement of the result. We shall always assume that k ^ 0: indeed, when k = we have 
Lq ~ Li — L2 and the spectral decomposition of this operator claimed in the latter lemma reduces to 
the Hankel transform. Since Lq = £2, it is obviously sufficient to consider the cases q = and q = 1. 

We shall define a new variable t := sin^ "^^j which ranges over (0, 1) if k > and over (—00, 0) if 
K < 0. Thus we can identify the Hilbert space ^^((O, diam(M)), d^s) with ^^((o, l), ^dt) (if k > 0) 
or L^((— 00, 0), ^dt) (if k < 0) via a unitary transformation V^, and write Lg = k V~^TqVn with 

To:=-t(l-<)^-(l-2i)^, Ti:=ro + ^. (33) 

It is not difficult to see that the singular differential operator Tq (resp. Ti) is in the limit circle 
case at and 1 (resp. at 0), and in the limit point case at infinity (resp. at 1 and at infinity). For 
simplicity of notation we shall still denote by Tg the self-adjoint operators defined by the action of the 
differential operators ([33|) on the domains [9] 

Dom(To) -.^ |w e iJ^((0, 1)) : Tqu e ^^((0,1)), \imtu'{t) ^ \im{l - t) u' (t) = o| , (34a) 

Dom(Ti) := |w G iJ^((0, 1)) : Tiu € L^((0,1)), \imtu{t) = o| (34b) 

when K > and 

Dom(ro) |w e iji((-cx),0)) : Tow e L^((-oo,0)), hmtu'(t) = o|, (34c) 

Dom(ri) := (u e iJi((-oo,0)) : Tiu e L2((-c»,0)), \imtu'(t)^o] (34d) 

when K < 0. The domains of the operators Lq defined in Section|4]are simply Dom(Lg) — V^^^(Dom(T^)). 

When S is compact, the spectrum of Lq is discrete and its spectral decomposition can be written 
as follows. 
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Lemma 7.1. Let us suppose that k > 0. Then Proposition 



holds with 



dpo(A)= ^<5(A-Kj(j + l))jdA, dpi(A)= ^<5(A-A^(/ + 3j + 2)) )dA, 

1/2 



Wo{s,Kj{j + 1)) = 



«;(2j + 1) 



47r 



Pj (cos -/ks) , 



wi{s,K{f + 3j + 2)) = 



<2j + 3) 



An 



1/2 



cos 



pf''^(cosV^s), 



where P^, and pjf'^^ respectively denote the Legendre and Jacobi polynomials of degree v . 

Proof. It is standard that the eigenvalues of Tq are nonnegative and that the normahzed eigenfunctions 
of Tq provide an orthonormal basis of L^{{0, 1)). The choice of boundary conditions (p4| ensures that 
the eigenvalues of Lq are the numbers A ^ for which the classical hypergeometric equation 

= (A - Tq)u = t{l - t)u" + (1 - 2t)u' + Au 

has a polynomial solution. The solution to the latter equation which is continuous at is proportional 
to the hypergeometric function ^ (i + -^Y^'^i 5 + (l ^" ^V^^^ 1; t) a-nd this fimction becomes a 



polynomial if and only if (4 + A) 



1/2 



Hence spec(ro) = + 1) : j € N}. In this case the 



latter hypergeometric function reduces to a Legendre polynomial, and using the well known formulas 
for the norm of these polynomials one readily arrives at the expression 

\/2jTTP,(l~2t) 

for the normalized eigenfunctions of Tg. By making use of the identification of L^{{0, diam(E)), d/is) 
with i^((0, 1), ^dt) via we immediately arrive at the above expression for the spectral resolution 

of Lq. 

The proof for Li is analogous. If we set u{t) =: (1 — t)v{t), the eigenvalue equation reads as 

= {A-Ti)u= (l-t) {t{l - t)v" + (1 - U)v' + (A - 2)v) , u G Dom(ri) . (35) 

The eigenvalues of Ti can be easily seen to coincide with the polynomial solutions to the latter equation. 
As the solution v of ([35]) regular at is proportional to P(| - {\+ A)^/^, | + (^ + A)^/^, 1; t), this 
impHes that {\ + A)^/^ - | G N and spec(Ti) = {p + 3 j + 2 : j e N}. Hence we immediately obtain 
the desired formula by writing the resulting eigenfunction function in terms of a Jacobi polynomial, 
expressing the result in the variable s and normalizing it to have unit norm. □ 



When S is noncompact, the spectrum of Lq is absolutely continuous and its spectral resolution 
can be derived using the Weyl-Kodaira theorem. 

Lemma 7.2. Let us suppose that k < and set 
Then Proposition \4-2\ holds with 



A. 



dpo(A) = dpi (A) = — tanh 
An 



1/2- 



dA, 



wo{s, A) — cosh 



-2a{\/K) 



F 



wi{s, A) 



cosh ■ 



-2a{\/K) 



F 



{a{X/K),a{\/K), 1; - 
{a{X/K) + 1,q;(A/k) 



tanh 



1. 1: -tanh^ 
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Proof. It is a standard spectral-theoretic result that the spectrum of Tq, which is is absolutely con- 
tinuous, is given by (—00, —j]- Hence let us take a complex number A with Re A ^ — ^ and nonzero 
imaginary part and consider the equation 

= (A - To)u = t{l - t)u" + (1 - 2t)u' + Ku , (36) 



where t £ R and Tq is to be interpreted as a formal differential operator. Then the solutions of 
satisfying the boundary condition required in Dom(To) at are proportional to 

w{t, A) := (1 - t)-"(^)F(a(A), a(A), 1; ^) , 

and those which are square integrable at infinity are proportional to 

w~{t,A) := (l-t)-"(^)F(a(A),a(A),2a(A); ^) 



if ImA < 0. 



w+(<,A) := (l-<)"(^)~iF(l-a(A),l-a(A),2-2a(A);^^) if ImA>0. 



For notational simplicity we shall henceforth write a instead of a (A). It should be noticed that the 
determination of the square root ensures that these functions depend analytically on A in the region 
Re A < — J. It is well known [1] that one can write e.g. as a linear combination of w and w~ as 
w+{t,A) = kiiA)w{t,A) + k2{A) w;-(i,A), with 

k,iA):=4^, k^iA):^ r(.mi - 2«) 



r(2a-l) ' ' r(l-a)2r(2a-l) 

Let us use the notation 

Wiu,v) := t{l - t) {u'{t)v{t) - u{t)v'{t)) 

for the reduced Wronskian of two solutions u and v of Eq. which is actually constant. The 

Weyl-Kodaira theorem can be used to prove that the self-adjoint operator Tq admits a spectral 
decomposition analogous that of Proposition 14.21 where w plays the role of the function Wq and dpq 
must be replaced by the absolutely continuous Borel measure 

dp(A) := ^ ^^(^^ dA 

2t:iW{w{-,A),w+{-,A)) 

on (—00, —\]- By T, 15.3.10] one has that the function A) can be written as 

w^{t,A) = -£i^^log(-t) + ^{t,A) , 
1(1 — ay 

where (p(-, A) is of class in a neighborhood of 0. As W{w{-, A), w+(-, A)) is constant, we immediately 
arrive at the formula 

H^(«.(,A),u;+(,A))-^(^"^") 



r(l-a)2 

for the reduced Wronskian, which in turn yields the expression 



dp(A) = tanh 



dA (37) 



for the measure p. If now use the identification of £^((0, 00), d/is) with L2((— oo,0), ^dt) we readily 
arrive at the formulas for the spectral decomposition of io- 

The analysis of Li is similar. Setting u{t) =: (1 — t)v{t) in the equation (A — Ti)w = as in 
Lemma [73] we find that the regular solution of this equation at is 



w{t,A) := (1 -t)-"F(a+l,a- 1, 1, ^-^j 
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up to a multiplicative constant, and that the solution which is square integrable at infinity is propor- 
tional to 



A) := (1 - ty^F{a + l,a- I, 2a 



1 



1 - t 



w+{t, A) := (1 - t)"-'^F(2 - a, -a, 2 - 2a; ^ 
Since the reduced Wronskian of w and is 

W{w{;A),W+i;A)) 

and 



- t 



if ImA < 0, 
if Im A > . 



r(2-2a) 



w+{t, A) 



r(Q + i)r(Q-i) 

r(2a- 1) 



w{t,A) 



T{2-a)T{-a) 
r(l- 2a) 



r(-a)r(2-a) 

the same reasoning as above shows that the spectral measure dp associated with Ti is also given by 
Eq. dSZl). □ 



As stated in Proposition 14.21 the integral representation (fH|) converges to g{Lq)u in the norm 
topology, and by Egorov's theorem this implies that the latter integral converges uniformly except 
on a subset of arbitrarily small measure. We find it convenient to conclude this section with another 
simple but useful observation concerning the pointwise convergence of previously defined the spectral 
decompositions. 

Lemma 7.3. Let Uq, Wq, Cq^^ o,nd pq he defined as in Proposition \4. 2\ and let u G C^([0, diam(E))). 
Then 



u{s) 



Wq{s,\) UqU{\) Apq{\) 



(38) 



pointwise for all s £ [0, diam(E)), and the convergence is uniform. 



Proof. For the sake of concreteness we restrict ourselves to the case k > 0. As C^([0, 1)) C 
Dom(L''") = V^r^(Dom(r™)) for any nonnegative integer m, it follows that the integral 



(L"u,«) 



\UqUi\)\^ dpq{X) < 



OO 



(39) 



must convergent for all m and u € C^([0, diam(S))). 

It is well known [I] that the Legendre and Jacobi polynomials satisfy 



p. 



(0,2) 



(01 ^ 



(J + 2)(J- + 1) 



for all ^ G [—1, 1]. By the explicit expression for Wq, this immediately implies that 



Wq{s,X)UqU{X) dpq{X) 



{ciX + C2y^^\UqU{X)\dpq{X) < 



which converges by virtue of Eq. ([39| . Here ci and C2 are positive constants. Hence the integral ([38]) 
converges uniformly if the support of u is contained in [0, if] with some K < diam(I]). (The same 
argument also yields convergence in the C'' strong topology) . When k ^ the proof is analogous and 
will be omitted. □ 
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